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We discuss production of charm quarks and mesons as well as nonpho- 
tonic electrons in pp scattering at RHIC. The distributions in rapidity and 
transverse momentum of charm and bottom quarks/ antiquarks are calcu- 
lated in the /cj-factorization approach. The hadronization of heavy quarks 
is done by means of fenomenological fragmentation functions. The semilep- 
tonic decay functions found by fitting semileptonic decay data measured by 
the CLEO and BABAR collaborations are used. Good description of the 
inclusive data at large transverse momenta of electrons is obtained and a 
missing strength at small transverse momenta of electrons is found. 

Furthermore we discuss kinematical correlations between charged lep- 
tons from semileptonic decays of open charm/bottom, leptons produced in 
the Drell-Yan mechanism as well as some other mechanisms not included so 
far in the literature. When calculating the Drell-Yan processes we include 
transverse momenta of q and q, using the Kwiecihski parton distributions. 
Reactions initiated by purely QED 7*7*-fusion in elastic and inelastic pp 
collisions as well as diffractive mechanism of exclusive cc production are 
included. A good description of the dilepton invariant mass spectrum 
measured by the PHENIX collaboration is achieved. Predictions for the 
dilepton pair transverse momentum distribution as well as distribution in 
azimuthal angle between electron and positron are presented. 

PACS numbers: 12.38.-t,12.38.Cy,14.65.Dw 

1. Introduction 

Some time ago the PHENIX and STAR collaborations have measured 
transverse momentum distribution of so-called nonphotonic electrons [H [2] • 
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It is a common wisdom that the dominant contribution to the nonphotonic 
electrons/positrons comes from the semileptonic decays of charm and/or 
beauty mesons. Formally such processes can be divided into three sub- 
sequent stages. First cc or bb quarks are produced. The dominant mecha- 
nisms being gluon-gluon fusion and quark-antiquark annihilation close to the 
threshold. Next the heavy quarks/antiquarks fragment into heavy charmed 
mesons D,D* or B,B*. The vector D* and B* mesons decay strongly pro- 
ducing (pseudo)scalar D and B mesons. Finally the heavy pseudoscalar 
mesons decay semileptonically producing electrons/positrons. 

The hadronization of heavy quarks is usually done with the help of 
phenomenological fragmentation functions, with parameters adjusted to 
the production of heavy mesons in e~^e~ or pp collisions. 

The last ingredient are semileptonic decays of heavy mesons. Only re- 
cently the CLEO [3] and BABAR [3] collaborations has measured precise 
spectra of electrons/positrons coming from the decays of D and B mesons. 
This is done by producing specific resonances: ^'(3770) which decays into 
D and D mesons (CLEO) and T(4S') which decays into B and B mesons 
(BABAR). In both cases the heavy mesons are almost at rest, so in practice 
one measures the meson rest frame distributions of electrons/positrons. 

In this presentation the results have been obtained within the fc^-factorization 
approach. At rather low RHIC energies intermediate x-values become rel- 
evant. The Kwiecinski unintegrated gluon (parton) distributions seem the 
best suited in this context [S]. We shall use both Peterson [9J and so- 
called perturbative [lOj fragmentation functions. The electron/positron de- 
cay functions fitted recently |12j to the recent CLEO and BABAR data will 
be used. 

Recently the PHENIX collaboration has measured dilepton invariant 
mass spectrum from to 8 GeV in pp collisions at ^/s = 200 GeV [7j. 
Up to now, production of open charm and bottom was studied only in 
inclusive measurements of charmed mesons [5J and electrons [B] and only 
inclusive observables were calculated in pQCD approach [TTl [T2]. Such 
predictions give rather good description of the experimental data, however, 
the theoretical uncertainties are quite large. 

Some time ago we have studied kinematical correlations of cc quarks [TJ] , 
which is, however, difficult to study experimentally. High luminosity and in 
a consequence better statistics at present colliders gives a new possibility to 
study not only inclusive distributions but also correlations between outgoing 
particles. Kinematical correlations constitute an alternative method to pin 
down the cross section for charm and bottom production. 

Below I shall limit to presentation of results obtained in |12[ [T5J . The 
original presentation at the conference included also diffractive processes. 
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2. Formalism 

Let us consider the reaction /ii + /i2 Q + Q + X , where Q and Q are 
heavy quark and heavy antiquark, respectively. 




Fig. 1. Basic diagrams relevant for gluon-gluon fusion (left panel) and quark- 
antiquark annihilation (right panel) in our fct-factorization approach. 



In the fct-factorization approach the multiply differential cross section 
reads: 



6^ {tii^t + fi2,t - Pu - P2,t) -Fi(a:i,Ki J J'j(x2,K|j , (1) 



where Fi{xi, k{ j) and Fj{x2, n\ t) are the so-called unintegrated gluon (par- 
ton) distributions. Leading-order matrix elements for off-shell gluons [15\ 
\T6\ [T7] were used. The two-dimensional Dirac delta function assures mo- 
mentum conservation. The unintegrated parton distributions are evaluated 
at: 

xi = ^ exp(yi) + ^ exp(y2), X2 = ^ exp(-?/i) + ^ exp(-?/2), where 

^i,t = \Jplt + ^Q- 

Introducing new variables: Qt = + K2,t , 
Qt = Kit — K2t can write: 



dai 



dyidy2d'^pi^td'^P2,t J 47r2 167r2(xiX2s)2 

^i{xi,i4,t) ^j{x2,i4^t) ■ (2) 

This formula is very useful to study correlations between the produced heavy 
quark Q and heavy antiquark Q |14) . 
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At the Tevatron and LHC energies the contribution of the gg — t- QQ 
subrocess is more than one order of magnitude larger than its counterpart 
for the qq — t- QQ subprocess. At RHIC energy the relative contribution 
of the qq annihilation is somewhat bigger. Therefore in the following we 
shall take into account not only gluon-gluon fusion process but also the 
quark-antiquark annihilation mechanism. 

The production of electrons/positrons is a multi-step process. The whole 
procedure of electron/positron production can be written in the following 
schematic way: 

TIT = TT72~ ® Dq^d ® fo^e , (3) 
dyd'^p dyd'^p 

where the symbol ® denotes a generic convolution. The first term respon- 
sible for production of heavy quarks/antiquarks. Next step is the process 
of formation of heavy mesons. We follow a phenomenological approach and 
take e.g. Peterson [9] and Braaten et al. [10] fragmentation functions with 
parameters from the literature [18] . The electron decay function accounts 
for the proper branching fractions. 

The inclusive distributions of hadrons can be obtained through a con- 
volution of inclusive distributions of heavy quarks/antiquarks and Q — )• h 
fragmentation functions: 

da-{yi,Pit,y2,P2t,ct>) _ f Dq^h{zi) Dq^h{z2) da{yi,p%y2,P2t,4>)^^^^^^ 



dyidp^tdy2dpgd(l) J zi Z2 dyidpf^dy2dp'^^u^ 

(4) 

where: = — , = — , where meson longitudinal fractions zi,Z2 G 
(0,1). 

We use decay functions fitted recently [12] to the CLEO and BABAR 
data. In our approach the electrons (positrons) are generated isotropically 
in the heavy meson rest frame. 



3. Results 

3.1. Single electron spectra 

Before we start presenting our results for the spectra let us focus for a 
moment on the decay functions discussed shortly above. In Figl2]we show 
our fit [12j to the CLEO and BABAR data. The good quality fit of the 
data will allow us to obtain reliable predictions for electron/positron single 
particle spectra. 

Now we shall concentrate on transverse momentum distribution of elec- 
trons/positrons measured recently by the STAR and PHENIX collabora- 
tions at RHIC [H [2]. In FiglHl as an example we show results obtained 
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p. (GeV) p, (GeV) 



Fig. 2. Our fit to the CLEO [3] and BABAR ^ data. 

with the Kwieciiiski UPDFs In Ref.|12j we have discussed in addition 
other UGDFs. Uncertainties due to different combinations of factorization 
and renormahzation scales as well as due to different choices fragmentation 
functions are shown in FigH] In these calculations we have included both 
gluon-gluon fusion as well as quark-antiquark annihilation. In the last case 
we use matrix elements with on-shell formula but for off-shell kinematics 
(the discussion of this point can be found in our earlier paper [14]). In 
Ref.[T2] we have discussed also uncertainties due to the choice of quark 
masses. 

Study of nonphotonic and hadron correlations allows to "extract" 
a fractional contribution of the bottom mesons B/{D -\- B) as a. function 
of electron/positron transverse momentum [19j. Recently the STAR col- 
laboration has extended the measurement of the relative B contribution to 
electron/positron transverse momenta ~ 10 GeV [20]. In FigjS] (Kwiecinski 
UPDFs) shown are results for different scales and different fragmentation 
functions. There is a strong dependence on the factorization and renor- 
mahzation scale in the case of the Kwieciiiski unintegrated gluon/parton 
distributions. A slightly better agreement is obtained with the Peterson 
fragmentation functions. 

3.2. Electron- positron correlations 

When calculating correlation observables we have included also processes 
shown in Figi6] and FiglT) The photon-photon induced processes were first 
included in Ref. |13j . The central exclusive diffractive process shown in 
FigJ7]was first proposed in Ref.[22j. 

In Figl5]we show e'^e~ invariant mass distributions calculated with the 
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Fig. 3. Transverse momentum distribution of electrons/positrons obtained with 
the Kwiecihski UPDFs. Different combinations of factorization and renormahza- 
tion scales are used. We show separately contributions of the gluon-gluon fusion 
(black) and quark- ant iquark annihilation (grey). On the left side results with the 
Peterson fragmentation functions and on the right side with BCFY fragmentation 
functions. 

Kwieclnski (left) and KMR (right) UGDFs. One can clearly see that both 
the Kwieclnski and KMR [21 J UGDFs give fairly good description of the 
data for Mg+g- > 3 GeV. At small invariant masses the Kwieclnski UGDF 
underestimates the PHENIX data and the KMR UGDF starts to overesti- 
mate the data points below Mg+g- = 2 GeV. 

In Figj9]we show uncertainties related to the contribution of semileptonic 
decays. Complementary the left panel presents uncertainties due to the 
factorization scale variation as described in the figure caption. The right 
panel shows uncertainties due to the modification of the heavy quark masses. 

If the detector can measure both transverse momenta of electron/positron 
and their directions, one can construct a distribution in transverse momen- 
tum of the dielectron pair: pt^sum = Pit +P2t- Our predictions including the 
semileptonic decays and Drell-Yan processes are shown in the left panel of 
FigJlOl Both processes give rather similar distributions. To our knowledge 
the distributions of this type were never measured experimentally. The dis- 
tribution in pt,sum is not only a consequence of gluon transverse momenta 
but invlolves also fragmentation process and semileptonic decays. With 
good azimuthal granulation of detectors one could also construct distri- 
bution in azimuthal angle between electron and positron. Corresponding 
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Fig. 4. Factorization and renormalization uncertainty band of our fcj-factorization 
calculation with unintegrated Kwiecinski gluon, quark and antiquark distributions 
for the Peterson fragmentation function (left panel) and BCFY fragmentation func- 
tion (right panel). The open triangles represent the PHENIX collaboration data 
and the solid circles the STAR collaboration data. 




p. (GeV) (GeV) 



Fig. 5. The fraction of the B decays for the Kwiecinski UPDFs. The uncertainty 
band due to the choice of the scales is shown for the Peterson (left) and Braaten 
et al. (right) fragmentation functions. Both gluon-gluon fusion as well as quark- 
antiquark annihilation are included in this calculation. 
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Fig. 6. Diagrammatic representation of processes initiated by photon-photon sub- 
processes: double-elastic, double-inelastic, inelastic-elastic and elastic-inelastic. 




Fig. 7. The mechanism of exclusive double-diffractive production of open charm. 




Mg.g (GeV) Mg.g (GeV) 

Fig. 8. Dielectron invariant mass distribution for pp collisions at \/s = 200 GeV for 
the Kwiecihski (left) and KMR (right) UGDFs. Different contributions are shown 
separately: semileptonic decay of charm by the blue solid line, semileptonic decay of 
bottom by the red solid line, Drell-Yan mechanism by the long dashed line, gamma- 
gamma processes by the blue dashed line and the central diffractive contribution 
by the green dotted line. In this calculation we have included azimuthal angle 
acceptance of the PHENIX detector [7]. 

predictions are shown in the right panel of FigllOi One can see an interest- 
ing dependence on the invariant mass of the dielectron pair - the smaller 
the invariant mass the large the decorrelation in azimuthal angle. 
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Mg.g (GeV) Mg-g (GeV) 

Fig. 9. The uncertainties of theoretical calculations. The left panel shows the 
factorization scale uncertainties, the lower curve corresponds to fJ.j?,li% — rn\ ^ + 
m\ f and the upper curve to = fcf , iJ^p = ^roq, where kt is gluon transverse 
momentum. The right panel shows the quark mass uncertainties as indicated in 
the figure. 




Fig. 10. Distribution in transverse momentum of the dielectron pair (left) and 
in azimuthal angle between electron and positron (right) for semileptonic decays 
(solid line) and Drell-Yan processes (dashed line). Here Kwiecinski UGDF and 
Peterson fragmentation function were used. 

4. Conclusions 

We have calculated inclusive spectra of nonphotonic electrons/positrons 
for RHIC energy in the framework of the fcj-factorization. We have con- 
centrated on the dominant gluon-gluon fusion mechanism and used unin- 
tegrated gluon distribution functions from the literature. Special emphasis 
was devoted to the Kwiecinski unintegrated gluon (parton) distributions. In 
this formalism, using unintegrated quark and antiquark distributions, one 
can calculate in addition the quark-antiquark annihilation process including 
transverse momenta of initial quarks/antiquarks. 

When calculating spectra of charmed (D, D*) and bottom (S, B*) 
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mesons we have used Peterson and Braaten et al. fragmentation functions. 
We have used recent fits to the CLEO and BABAR collaborations data for 
decay functions of heavy mesons. 

Our results have been confronted with experimental data measured re- 
cently by the PHENIX and STAR collaborations at RHIC. A reasonable 
description of the data at large transverse momenta of electrons/positrons 
have been achieved. We have discussed uncertainties related to the choice 
of the factorization and renormalization scales and those related to the frag- 
mentation process. Although the uncertainty bands are rather large, there 
is a missing strength at lower electron/positron transverse momenta. 

We have discussed also correlations of charmed mesons and dielectrons 
at the energy of recent RHIC experiments. We have calculated the spec- 
tra in dielectron invariant mass, in azimuthal angle between electron and 
positron as well as for the distribution in transverse momentum of the pair. 
The uncertainties due to the choice of UGDFs, choice of the factorization 
and renormalization scales, have been discussed. We have obtained good 
description of the dielectron invariant mass distribution measured recently 
by the PHENIX collaboration at RHIC. 

At RHIC the contribution of electrons from Drell-Yan processes is only 
slightly smaller than that from the semileptonic decays. The distributions 
in azimuthal angle between electron and positron and in the transverse 
momentum of the dielectron pair from both processes are rather similar. 
It was found that the distribution in azimuthal angle strongly depends on 
dielectron invariant mass. 

We have also included exclusive central-diffractive contribution discussed 
recently in the literature. At the rather low RHIC energy it gives, however, 
a very small contribution to the cross section and can be safely ignored. 

The QED double-elastic, double-inelastic, elastic-inelastic and inelastic- 
elastic processes give individually rather small contribution but when added 
together are not negligible especially at low dielectron invariant masses 
where some strength is missing. 
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